Phase Structure in a Hadronic Chiral Model 
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We study the phase diagram of a hadronic chiral fiavor-SU(3) model. Heavy baryon resonances 
can induce a phase structure that matches current results from lattice-QCD calculations at finite 
temperature and baryon density. Furthermore, we determine trajectories of constant entropy per 
net baryon in the phase diagram. 
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Understanding hot and baryon-dense QCD matter is 
of central importance in theoretical and experimental 
heavy-ion physics. Various effective theories for chiral 
symmetry restoration predict that a line of first-order 
phase transitions in the plane of quark-chemical poten- 
tial fj, q versus temperature T ends in a critical point as 
the chemical potential is lowered (see 1] for a review). 
Presently, the lattice locates that point at T w 160 MeV, 
fi q w 120 MeV 0. The Compressed Baryonic Matter 
(CBM) experiment at GSI FAIR is planned to perform a 
dedicated experimental effort to detect that line of first- 
order phase transitions in relativistic heavy-ion collisions. 
It is hoped that by varying experimental parameters like 
the beam energy one could trigger phase transitions of 
variable strength (latent heat) and perhaps even locate 
the expected second-order critical point. 

Models relying exclusively on order-parameter dynam- 
ics typically predict significantly lower chiral phase tran- 
sition temperatures in baryon-dense matter than those 
found on the lattice (see e.g. Fig. 6 in 111). As shown 
by Gerber and Leutwyler some time ago |3(], while heavy 
hadronic states are suppressed by the Boltzmann factor, 
their contribution to the energy density at high temper- 
ature is substantial. This agrees with other studies us- 
ing a hadron resonance gas approach, which provides a 
reasonable description of the thermodynamics obtained 
on the lattice below the critical temperature 0. Heavy 
states also reduce the strong dependence of the "crit- 
ical temperature" (defined via the peak of a suitable sus- 
ceptibility) on the pion mass obtained in simple models 
for chiral order parameter dynamics [jj. The lattice in- 
dicates a relatively weak dependence of T c on the pion 
mass |6J. 

In this Letter we investigate the role of heavy hadronic 
states on the location of the chiral critical point within a 
non- linear SU(3) L x SU{3) R chiral model 0. Here, the 
phase transition at high temperature and baryon density 
is "driven" by baryonic resonance degrees of freedom, 
as suggested by the discussion above. We shall show 
that the model is able to reproduce not only the sketched 
qualitative phase structure but also the location of the 
endpoint. The properties of the high mass states (masses 
and couplings) are important for the actual location of 



the chiral phase transition line || in the plane of T and 
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Lattice results show that the susceptibility peaks of 
the chiral condensate and of the Polyakov loop coincide 
at /j, q = |{J which indicates that for small fi q those 
transition(s) involve a coupling of the chiral dynamics to 
the gauge fields, see e.g. 0]. However, within a ma- 
trix model for Polyakov loops, the effect of \i q > on 
the critical temperature for deconfinement is suppressed 
by 1/N C [ll|. I n contrast, /^-effects on the chiral critical 
temperature could be relatively strong. Indeed, models 
which couple chiral fields to Polyakov loops do suggest 
a decoupling of chiral and deconfinement dynamics at 
large chemical potential, and that the largest increase in 
the Polyakov loop still appears at T w T c (fi = 0), while 
the peak of the chiral susceptibility is shifted to signifi- 
cantly smaller T Thus at large chemical potential 
a confined phase with partially restored chiral symmetry 
might exist. 

Our chiral hadronic SU (3) Lagrangian incorporates 
the complete set of baryons from the lowest flavor-S'C/ (3) 
octet, as well as the entire multiplets of scalar, pseu- 
doscalar, vector and axialvector mesons In mean- 

field approximation, the expectation values of the scalar 
fields relevant for symmetric nuclear matter correspond 
to the non-strange and strange chiral quark condensates, 
namely the a and its ss counterpart £, respectively, and 
further the ui and <fi vector meson fields. Another scalar 
isoscalar field, the dilaton x, is introduced to model the 
QCD scale anomaly. However, if \ does not couple 
strongly to baryonic degrees of freedom then it remains 
essentially "frozen" below the chiral transition. Conse- 
quently, we focus here on the role of the quark conden- 
sates. 

Interactions between baryons and scalar (BM) or vec- 
tor (BV) mesons, respectively, are introduced as 

4m = - ^2 A (9ta<y + g^C) ipi , (1) 

i 

£bv = -y^t (gtujjou + ff#7o0°) ipt , (2) 

i 

Here, i sums over the baryon octet (TV, A, E, S). A term 
£ VO c with mass terms and quartic self-interaction of the 



2 



vector mesons is also added: 
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The scalar self-interactions are 
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Interactions between the scalar mesons induce the spon- 
taneous breaking of chiral symmetry (first line) and the 
scale breaking via the dilaton field \ (\&st two terms). 

Non-zero current quark masses break chiral symmetry 
explicitly in QCD. In the effective Lagrangian this corre- 
sponds to terms such as 
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According to £bm DJ, the effective masses of the 
baryons, m*(cr, C) = <JW cr + g^ £ , are generated through 
their coupling to the chiral condensates, which attain 
non-zero vacuum expectation values due to their self- 
interactions in £ ©■ The effective masses of the 
mesons are obtained as the second derivatives of the 
mesonic potential about its minimum. 

The baryon-vector couplings g^ and result from 
pure /-type coupling as discussed in 0, g^ = (n q — 

n \)9& > 9i4> = ~( n l " n\)y/2gX , where g% denotes the 
vector coupling of the baryon octet and n l the num- 
ber of constituent quarks of species i in a given hadron. 
The resulting relative couplings agree with additive quark 
model constraints. 

All parameters of the model discussed so far are fixed 
by either symmetry relations, hadronic vacuum observ- 
ables or nuclear matter saturation properties (for details 
see 0). In addition, the model also provides a satisfac- 
tory description of realistic (finite-size and isospin asym- 
metric) nuclei and of neutron stars El • 

If the baryonic degrees of freedom are restricted to the 
members of the lowest lying octet, the model exhibits a 
smooth decrease of the chiral condensates (crossover) for 
both high T and high /j, 00] . However, additional bary- 
onic degrees of freedom may change this into a first-order 
phase transition in certain re gime s of the T-fi q plane, de- 
pending on the couplings [M TL2I Ilij . To model the in- 
fluence of such heavy baryonic states, we add a single 
resonance with mass tub, — too + gna and vector cou- 
pling gnu, — rygNuj ■ The mass parameters, too , gR and 
the relative vector coupling ry represent free parame- 
ters, adjusted to reproduce the phase diagram discussed 
above |l9j. In principle, of course, one should include 
the entire spectrum of resonances. However, to keep the 
number of additional couplings small, we effectively de- 
scribe all higher baryonic resonances by a single state 
with adjustable couplings, mass and degeneracy. 



In what follows, the meson fields are replaced by their 
(classical) expectation values, which corresponds to ne- 
glecting quantum and thermal fluctuations. Fermions, of 
course, have to be integrated out to one-loop. The grand 
canonical potential can then be written as 



fl/V = —C VC c ~ Co ~ £sb — V v . 
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where 73,771/ denote the baryonic and mesonic spin- 
isospin degeneracy factors, respectively, and E* B M (k) = 

k 2 + m* B M 2 are the corresponding single particle en- 
ergies. In the second line we sum over the baryon 
octet states plus the additional heavy resonance with de- 
generacy n (assumed to be 16). The effective baryon- 
chemical potentials are /J,* = fa — g^uj — g^ifi, with 
Hi = (n q — n l q )n q + (n l s — n\)n s . The chemical poten- 
tials of the mesons are given by the sum of the corre- 
sponding quark and anti-quark chemical potentials. The 
vacuum energy V vac (the potential at ps = T — 0) has 
been subtracted. 

By extremizing Q/V one obtains self-consistent gap 
equations for the mesons. Here, we consider non-strange 
matter, adjusting /i s for any given T and \i q so that the 
net number of strange quarks is zero. The dominant 
"condensates" then are the a and the u> fields. 

The properties of the QCD transition at vanishing 
chemical potential depend on the number of quark flavors 
and on their masses. In quenched lattice gauge theory, 
which corresponds to infinitely heavy quarks, the criti- 
cal temperature is T c w 260 MeV |l5j. Smaller quark 
or pion masses, respectively, reduce the critical tempera- 
ture . Lattice calculations with 2 and 3 flavors of light 
quarks predict critical temperatures of T c « 175 MeV 
and T c « 155 MeV, respectively. However, we note 
that recently 0] values as high as 190 MeV were re- 
ported for the 2-flavor case with reasonably small quark 
masses, and for reduced pion over-degeneracy on the 
lattice. For two flavors, the lattice results support a 
continuous transition in agreement with symmetry ar- 
guments 0]. For three massless flavors, a first-order 
phase transition is expected, which is washed out at a 
pion mass of m w « 70 - 260 MeV flU ]. 

Within our model, at zero baryon density we obtain 
values for the critical temperature around 140— 190 MeV, 
which are in reasonable agreement with the numbers ob- 
tained from lattice calculations. Furthermore, we observe 
a rather slow increase of the critical temperature with 
increasing explicit mass term too (corresponding to the 
quark mass term of QCD, cf. Fig. GJ. Finally, the first- 
order phase transition observed for small explicit mass 
term indeed turns into a crossover beyond some criti- 
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FIG. 1: Critical temperature T c at vanishing chemical po- 
tential vs the explicit mass term mo for different values of 
the scalar coupling gn. Full lines denote a first-order phase 
transition and dotted lines correspond to a crossover. 

cal value of mo- Neglecting the octet hyperons, which 
basically corresponds to the two- flavor limit, leaves a 
crossover at fi = and increases the critical temperature. 
If in turn we assign the additional resonance to a spin-3/2 
decuplet of flavor-S'C/ (3) and take the 3-flavor limit then 
chiral restoration happens via a first-order phase transi- 
tion and the critical temperature drops by « 30 MeV. 
This behavior is also in (qualitative) agreement with the 
lattice QCD findings described above. 

We now turn to non-zero (net) baryon density. Within 
certain ranges for the mass and couplings of the heavy 
resonance, the location of the critical endpoint is reason- 
ably close to the lattice results. This is in contrast to 
models which do not account for heavy degrees of free- 
dom and usually predict the endpoint at too low temper- 
ature (see e.g. |l| for a summary). Choosing, for example, 
w_r(cto) = 2 GeV, ry — 0.4 and m = 0.57 GeV results 
in T E w 180 MeV |see Fig. 0). This is a little higher 
than Te from Ref. [2| but the peak of the susceptibility 
of the order parameter at /i = is at T c — 185.6 MeV, 
in good agreement with the recent lattice results at zero 
density. Smaller explicit mass term ~ mo for the addi- 
tional heavy resonance moves the critical endpoint closer 
to the temperature axis, finally turning the crossover at 
fi q = into a first-order transition (cf. Fig.QJ if mo drops 
below a critical value. We also point out that the phase 
transition at T = occurs at a density above the nuclear 
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FIG. 2: First-order phase transition line in the T-p plane. 
m R (a ) = 2 GeV, r v = 0.4 and m = 0.57 GeV. 

Figure [3] shows the density and entropy contributions 
to the latent heat along the transition line of Fig. [21 The 
latent heat first grows with increasing chemical poten- 
tial. However, around p q = 150 MeV it decreases again, 
exhibiting a local minimum at (j, q w 200 MeV before in- 



creasing again towards low temperature. This structure 
is also present in the entropy contribution. The density 
contribution to the latent heat up to p, q « 200 MeV is 
flat but then rises very rapidly with increasing chemical 
potential. 
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FIG. 3: Top: Latent heat and its entropy (TAs) and density 
(/iflAps) contributions along the phase transition line versus 
the quark chemical potential. Bottom: Difference of the scalar 
density (Apl) in the two phases at the phase transition versus 
the quark chemical potential. 

This somehow unexpected behaviour is due to a change 
in the dominant degrees of freedom. Around p q « 
100 — 150 MeV nucleons dominate the discontinuity of 
the scalar density (cf . Fig. 131 at the phase boundary, but 
this contribution decreases above \x q « 150 MeV and 
finally vanishes at fi q w 270 MeV. In contrast, the con- 
tribution of the resonances increases sharply around 200 
MeV. Both effects together generate the small dip in the 
latent heat and the 'nose' in the e-ps phase diagram, 
cf. Fig. |SJ This behaviour is actually not that surprising 
given that lattice results at p q w suggest just one single 
peak in all susceptibilities. In contrast, at zero tempera- 
ture the jump from zero to finite nucleoli density signals 
the liquid-gas phase transition while the chiral transition 
is expected to happen at larger chemical potential [l7| . 
driven by newly populated heavier degrees of freedom. 
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FIG. 4: Phase transition line and adiabats in the T-\i q phase 
diagram. 
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Figure ^ depicts lines of constant entropy per net 
baryon (adiabats) in the phase diagram. These corre- 
spond to perfect-fluid expansion trajectories. At the 
phase transition they "bend" slightly towards smaller 
chemical potential. For the present EoS, the S/A « 20 
adiabat goes right through the endpoint of first-order 
transitions, higher specific entropies then correspond to 
crossovers. To make contact with experiments, e.g. at 
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FIG. 5: Adiabats in the e-pB phase diagram. The grey area 
represents the mixed phase. 

the upcoming GSI-FAIR facility, we use a simple esti- 
mate for the initial baryon number and energy density: 
p l ™ = 2 7cmPo, e™ = -\/spo7CM- This assumes that 
the entire initial beam energy and baryon number equi- 
librates in a Lorentz-contracted volume determined by 
the overlap of projectile and target in the center-of-mass 
frame. In the energy range of interest here, this simple 
model reproduces the specific entropy from more realis- 
tic three- fluid models quite well 18] . Within this sim- 
ple model beam energies between 5-10^4 GeV are suffi- 
cient to overshoot the shaded phase coexistence region 
(Fig. [3J). In the energy range from 5-25 A GeV, adia- 
batic expansion leads to a first-order phase transition 
back to the symmetry broken state with a latent heat 
of « 120 — 160 MeV/fm 3 . Higher collision energies lead 
to weaker (more "critical") first order transitions. The 
critical end point of the line of first-order transitions is 
reached for a beam energy of about 60 AGeV, with higher 
energies leading to crossover transitions. 

In summary, we have explored the influence of heavy 
baryon resonances on the phase diagram of a flavor- 
SU(3) hadronic chiral Lagrangian. We find that, with 
appropriate couplings, the model is able to reproduce 



the structure of the QCD phase diagram at finite den- 
sity in that a line of first-order phase transitions at high 
density terminates in a critical point at Te ~ 180 MeV, 
/i£ « 110 MeV. The interplay of nucleon and resonance 
contributions to the density may lead to a non-monotonic 
latent heat along the phase boundary. 
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Note that instead of an explicit mass term we could have 
coupled the resonance to the dilaton \ . 



